Shock wave in Euler - Heisenberg - Kockel nonlinear vacuum. 
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An analytical approach to the theory of electromagnetic waves in nonlinear vacuum is developed. The 
evolution of the pulse is governed by a system of nonlinear wave vector equations. Exact solution with own 
angular momentum in form of a shock wave is obtained. © 2012 Optical Society of America 
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Contemporary hight-power laser facilities can gener- 
ate optical pulses with intensities of order 10 22 W/cm 2 . 
At the same time the critical power for observation self- 
action effects due to virtual electron-positron pairs is 
of order [MS] P cr = X 2 /8n n 2 = 2.5 - 4.4 x 10 24 W, 
at a wavelength 1 fim. Thus, for a laser pulse with 
waist r_i_ = 1 mm the corresponding intensity becomes 
/™ c = P cr /r\ = 2.5 - 4.4 x 10 26 W/cm 2 , which is above 
the range of the new high-power lasers. The nonlinear ad- 
dition to the refractive index in vacuum depends also on 
the magnetic field. That is why new different nonlinear 
effects can be expected. There are not only self-action ef- 
fects, but also vacuum birefringence [H6], different kinds 
of four wave interaction 0[7l[8] and higher order har- 
monic generation [9]. In this paper we shall investigate 
the self-action effect only for intensities of order I^ c . 

Euler, Heisenberg and Kockel PUM] predict one in- 
trinsic nonlinearity of the electromagnetic vacuum due 
to the electron-positron nonlinear polarization. The clas- 
sical field-dependent nonlinear vacuum dielectric tensor 
can be written in the form 
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where a complex form of presenting of the electrical Ei 
and magnetic Bi components is used. Note that the term 
containing BiB^ vanishes, when a localized electromag- 
netic wave with only one magnetic component Bi is in- 
vestigated. The dielectric response relevant to such op- 
tical pulse is thus 
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In the case when the spectral width of a pulse Afc 2 ex- 
ceeds the values of the main wave- vector, i.e. Ak z ~ ko, 
the system of amplitude equations can be reduced to 
wave type [T2] and in nonlinear vacuum becomes 
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where 7 



and E, B are the amplitude func- 
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tions. Initially, we can write the components of the elec- 
trical and magnetic fields as a vector sum of circular 
and linear components E z ; E c = iE x — E y ; B\ = —B z . 
Thus ((3]) is transformed in the following scalar system of 
equations 
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Let us now parameterize the 3D + 1 space-time 
through pscudospherical coordinates (r,T,8,ip): 

z = r cosh(r) cos(#), y = r cosh(r) sin(0) sin(</?), x = 
?*cosh(r) sin((9) cos(<ys) and ct — rsinh(r), where r = 
y/ 'x 2 + y 2 + z 2 — c 2 t 2 . After calculations the correspond- 
ing d'Alambert operator in pscudospherical coordinates 
becomes [13] 
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where with Ag iV is denoted the angular part of the usual 
Laplace operator 
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The system of equations Q in pseudo-spherical coordi- 
nates becomes 
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Fig. 1. Time evolution of the intensity profile / of the spherically symmetric analytical solution (|17[) of the linear 
wave equation (fT5|) (ro = 1 and c = 1). The initially (t = 0) localized amplitude function (Fig. la) decreases with 
the generation of outside and inside fronts (Fig. lb), while the energy density distributes over the whole space for a 
finite time (t = 10). 
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Eq's ([7]) are solved using the method of separation of the 
variables. 



Ei(r,T,6,<p) = R(r)Ti(T)Yi(6,<p) 
B l (r,T,0,<p)=R(r)T l (T)Y l (0,<p), 



(8) 



where i = z,c. We use an additional constrain on the 
angular and "spherical" time parts 



\T Z \ 2 \Y Z \ 2 



\T C \ 2 \Y C \ 2 - \Ti\ 2 \Yi 



2.VJI2 



const. 



(9) 



The condition ([5]) separates the variables. The nonlinear 
terms appear in the radial part only. Thus the radial 
parts obey the equation 



where a, 7 and the separation constants Ai, i = z,c,l 
satisfy the relations a 2 — 1 = Af, 2a 2 = 7. The cor- 
responding r - dependent part of the equations (J7J) arc 
linear 
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where i = z, c, Z and Ci are another separation constants 
connected with the angular part of the Laplace opera- 
tor Yi(6, if). Only the following solutions of the Eq. (fT2j) 
exist which satisfy the condition T z = coshr; T c — 
coshr; Ti = sinhr, with separation constants: for the 
electrical part A z = A c = 3; C z = C c = 2 and for the 
magnetic part A\ = 3; C; = 0. Thus the magnetic part 
of the system of equations ([7]) do not depends on the 
angular components, i.e. Yi{9,<f) = 0, as for the electri- 
cal part Y z {9,(f), Y c (9,tp) we have the following linear 
system of equations 
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where now i = z,c. There are only two solutions of the 
Eq. (|T3)) which satisfy the condition (J9j> : Y z = cos 8; Y c = 
sin6*exp(iiy9), Using the relation between the separation 
constants Ai and the real number a we obtain the Mow- 
ing values for a and 7: a 2 = 4; a = ±2; 7 = 8. 

Finally, we can write the exact solution of the system 
of nonlinear equations Q which describes the propaga- 
tion of a electromagnetic wave in nonlinear vacuum 
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where Ai, i = z,c, I are separation constants. We look 
for solutions which possess more clearly expressed local- 
ization than the scalar soliton solution obtained in [T2] 
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If we rewrite the solution in Cartesian coordinates, it is 
(11) not difficult to show that the solution ([M]) of the system 
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Fig. 2. Evolution of the intensity profile (jf 6|) of the solution (|15|) of the nonlinear vector wave equation in vacuum 
(c = 1) for i = and t = 10 correspondingly. The nonlinear wave demonstrates entirely different evolution, than the 
linear spherical ones: the shock wave preserves its amplitude maximum and self-compresses in r direction. 



(0| admits finite energy and the electrical part possesses 
angular momentum I = 1 
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where r = y x 2 + y 2 + z 2 — c 2 t 2 . The intensity profile of 
the solution now becomes 
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For a comparison, on Fig. la and Fig. lb is shown the 
time evolution of the intensity profile / of a spherically 
symmetric analytical solution 
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of the linear scalar wave equation 
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obtained recently by us applying the Fourier method. 
We have used normalized scales ro = 1, c = 1, times 
of evolution t = and t = 10. The initially localized 
amplitude function in the linear case decreases with the 
generation of outside and inside fronts, while the energy 
density distributes over the whole space for a finite time. 
The evolution of the intensity profile (|16[) is presented 
on Fig. 2a and Fig. 2b for t — and t = 10 correspond- 
ingly. It is clearly seen from Fig. 2 that the solution (fTB"j) 
describes a nonlinear shock wave in vacuum. The wave 
admits entirely different evolution than the linear spher- 
ical ones: as the linear wave front enlarges spherically, 
the shock wave preserves its amplitude maximum and 
self-compresses itself in r direction. 

In this paper the nonlinear vector wave equations in 
nonlinear vacuum (|2|) is solved through the method of 



separation of the variables in a pseudo-spherical coordi- 
nate system. The obtained analytical solution (fT5")l rep- 
resents a spherical shock wave with its own angular mo- 
mentum I ~ 1 for the electrical field. Such high in- 
tensity wave can be generated not only from the laser 
sources, but also in a nuclear reaction, where a non- 
linear polarization of virtual electron-positron pairs ap- 
pears at the beginning. If wc compare the nonlinear vac- 
uum shock wave with a spherically symmetric solution 
of the linear wave equation, the difference becomes obvi- 
ous. While the spherically symmetric solution of the lin- 
ear wave equation forms inside and outside wave fronts 
and the amplitude significantly decreases, the nonlinear 
shock wave preserves the amplitude maximum and sclf- 
comprcss in r direction. 
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